Theory and data analysis for excitations in liquid ^He beyond the roton minimum 
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The hybridization of the single-excitation branch with the two-excitation continuum in the mo- 
mentum region beyond the roton minimum is reconsidered by including the effect of the interference 
term between one and two excitations. Fits to the latest experimental data with our model allow 
us to extract with improved accuracy the high momentum end of the ^He dispersion relation. In 
contrast with previous results we find that the undamped excitations below two times the roton 
energy survive up to Q = 3.6 due to the attractive interaction between rotons. 
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Although excitations in superfluid "^He have been 
widely studied in the last decades (see for instance 
Ref. 10,^), the nature of the single particle spectrum ter- 
mination remains unclear. Forty years ago, Pitaevskii 
predicted different kinds of termination depending on the 
detailed form of the spectrum at low momenta . In the 
case of a decay into pairs of rotons, Pitaevskii theory de- 
scribes the avoided crossing of the bare one-excitation 
branch with the continuum of two excitations. Within 
this picture the low energy pole is repelled by the con- 
tinuum, so that the spectrum flattens out for large Q to- 
wards 2 A losing spectral weight (A is the roton energy). 
At the same time, a damped excitation for u > 2A ap- 
pears and shifts to higher energies. Neutron scattering 
experiments later suggested that the decay of excitations 
into pairs of rotons did actually take place for momen- 
tum Q > 2.6 The flattening of the spectrum at 
an energy of the order of 2A being the main feature in 
support of Pitaevskii's picture. Despite the good quali- 
tative agreement between theory and experiment, several 
issues addressed by the theory could not be verified due 
to insufficient instrumental resolution. In particular, the- 
ory predicts a singular termination of the spectrum at a 
definite value of momentum Qc for a repulsive interac- 
tion (V4 > 0) between rotons. In the case of an attrac- 
tive interaction, hybridization between the roton bound 
state and the single quasiparticles is instead expected, as 
proposed by Zawadowski-Ruvalds-Solana (ZRS) ||], with 
the consequence that an undamped quasiparticle peak at 
an energy slightly below 2 A should be present [||-§|. To 
our knowledge it has not yet been possible to distinguish 
clearly between these two cases. 

Smith et al. Q performed the first complete experi- 
mental analysis for 2.7 < Q < 3.3 A~^. Although they 
found indications of a repulsive interaction (V4 > 0) using 
ZRS theory, they pointed out that the experimental find- 
ing of energies for the quasiparticle peak above 2A was 
not accounted for by the theory with reasonable values of 
the parameters. As a matter of fact, the position of the 
low-energy peak was not extracted from the data with 
ZRS theory, but rather by fitting a Gaussian peak on 
a background of constant slope. The resulting spectrum 
reached values above 2A, in contrast with theoretical pre- 



dictions. From the theoretical point of view it is not pos- 
sible to explain in a simple way the presence of sharp 
peaks at energies above 2A, as the corresponding excita- 
tion should be unstable towards decay into two rotons. 
Moreover, the experimental finding of a repulsive interac- 
tion disagrees with different theoretical calculations that 
predict a negative value for V4 (l^J^. More recent ex- 
perimental investigations by Fak and coworkers 12 
concentrate on the interesting temperature dependence 
of the dynamical structure factor S'(Q,a') and do not ad- 
dress directly the issue of the quasiparticle energy or the 
interaction potential between rotons. In any case they 
show clearly that there is a strong correlation between 
the low-energy peak and the high energy continuum as 
Q increases from 2.3 to 3.6 A~^ |l^, thus indicating that 
hybridization takes place. We further recall that this 
hybridization is expected as direct consequence of Bose 
condensation as explained in details in Ref. 

From the theoretical point of view, we recall that the 
validity of Pitaevskii and ZRS theories is restricted to a 
small region around 2A. Indeed Pitaevskii in his orig- 
inal paper [|j exploited the logarithmic divergence ap- 
pearing in the two-roton response function [_F'o(Q,cl') ~ 
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Lu{Q), lli{Q) is the measured spectrum and p = (Q, lu)] to 
solve exactly the many-body equations. This elegant the- 
ory provides explicit expressions for the Green and the 
density-density correlation functions, valid only in the 
small energy range where the singularity dominates. This 
fact leads to problems in data analysis when Q grows so 
that the bare excitation energy uJoiQ) takes values above 
2A. In fact, the signal around 2A strongly decreases with 
Q and spectral weight shifts to higher energies following 
LUo{Q). To understand the correlation between the high 
energy part of the spectrum and the onc-excitation con- 
tribution, it is necessary to extend the validity of the 
theory to a wider range of energies in order to describe 
properly the continuum contribution to S{Q,uj). It then 
becomes crucial to consider the effect of the direct excita- 
tion of two quasiparticles by the neutron and its interplay 
with the one-quasiparticle excitations usually considered. 

The aim of the present paper is to construct such an 
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extension of the theory to describe experimental data for 
S'(Q,w) at very low temperatures. Since excitations in 
■^He are stable (r(Q)/tj(Q) « IQ-^ for rotons at 1.3 K 
with r the half width of the excitation), it is an excel- 
lent approximation to write the Hamiltonian directly in 
terms of the creation and destruction operators and b 
of these excitations: 
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We consider for the moment only the V3 interaction that 
induces the hybridization of the single with the double 
excitation. In general this vertex will be a function of 
two momenta, for instance the total momentum of the 
two particles and the momentum of one of them. We ne- 
glect from the outset this dependence on momentum as 
it is expected to be smooth in the region of interest and 
less important than the frequency dependence retained 
in the following. Since S{CI,lu) is the imaginary part of 
the density-density response function [x(Q,a;)] it is con- 
venient to express the density operator pp in terms of the 
b and 6^ field operators. In general this will be an infinite 
series in the ^-fields, by retaining only the one and two 
quasi-particle terms we obtain: 
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Eq. (H) gives the most general second order form for pp 
in terms of b and 6^ that fulfills parity, time-reversal, and 
p^p = p-p transformation properties. For the same in- 
variances the three functions introduced a, (3, and 7 are 
bound to be real. Although the above expression for pp 
is quite general it can be obtained microscopically within 
a particular approximation scheme [|l3| 
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FIG. 1. Full lines stands for the exact Green function, 
dashed boxes for the sum of all one-particle irreducible di- 
agrams, and the black box for the sum of all diagrams with 
two lines closing at the two ends. x(p) is the sum of all dia- 
grams with two density insertions (wiggly lines). 



It is possible at this point to calculate perturbatively 
x(Q,w) using the explicit form of pp given in Eq. (g) 
and the Hamiltonian (|l|). At zero temperature the con- 
tribution of the 7(p, q) term to x vanishes, while the 
momentum dependence of a(p) and /3(p,q) is neglected 
for the sake of simplicity. The diagrammatic theory for 
the model is shown in Fig. ^ together with the definition 
of x{p), Fo{p), Ftip), F{p), and Goip) = [oJ - Wo(Q)]-^ 
In this figure dashed boxes represent the sum of all one- 
particle irreducible diagrams. Ft stands for the sum of all 
diagrams with two lines closing at the two ends, and F is 
the self-energy without the two external couplings. Since 
the momentum dependence of the vertices is neglected, 
the Dyson equations for Q and Ft (see Fig. ||) take the 
simple form of a coupled algebraic system: 

Q{p) = Qo{p) + Qo{p) ^3 F{p) g{p), (3a) 
Ft{p) = F{p) + F{p) Goip) V3 Ft{p). (3b) 

By solving Eqs. (^) and substituting Q and Ft into the 
expression for x{p) we finally obtain: 
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Eq. (^ is the basic expression that we will use in the fol- 
lowing for the fits to the data [5(Q,a;) = — Imx(Q,a;)]. 
The presence of a V4 interaction does not change the 
above treatment since no V4 vertex can connect F to Q 
at zero temperature. This implies that the additional dia- 
grams due to a V4 interaction contribute only to F . Since 
the explicit calculation of F is difficult and in general de- 
pends on the detailed structure of the vertex functions, 
we prefer to extract it directly from data by exploiting 
the large (energy and momentum) region of validity of 
Eq. (^. We proceed by noting that in Eq. (||) only F 
and uJo (through Qo) depend on Q. Explicit evaluation of 
Fo suggests that the momentum dependence of F should 
not be pronounced in the range of interest 2.3-3.2 A"!. 
We thus drop this dependence completely in F and fit the 
resulting function F{lo) to the data. The Q-dependence 
of iOo is not negligible because it drives uJo(Q) through 
the 2A line. In this way we are left with only one param- 
eter dependent on Q. It is now possible to extract both 
uj{Q) and F{ijj) by fitting Eq. (^) to all sets of data of 
different momentum at the same time. This procedure 
exploits fully the information contained in the data be- 
cause it is sensitive to the correlation among sets with 
different Q. For this reason we are able to extract infor- 
mation for w(Q) and for the function F{uj) over a scale 
(slightly) smaller than the instrumental resolution. 

A few words are necessary to explain how we can find 
the complex function F{lo) from the data. The imag- 
inary and real part of F are related by the relation 
ReF(w) = -I/ttP / dLo'lTa.F{uj)/{u} - lo'). We thus need 
only to parametrize one of them, we choose Imi^(a;), since 
on this function it is easier to apply the physical con- 
straint that excitations are kinematically stable for ener- 
gies smaller than 2A. This condition reads ImF(ijj) = 
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for ui < 2 A (the small spectral density between a;(Q) and 
2 A can be neglected). A simple way to parametrize the 
function ImF is then to choose a set of values of ui say 
{wi, UJ2, ■ ■ ■ , ujn} with 2A — uji < uj2 < ■ ■ ■ < reason- 
ably spaced and to assign a free parameter, a^, to each 
of them. Imi^(a;) can then be defined as a cubic spline 
interpolation on such a set. The integral to calculate 
the real part can be performed analytically so that its 
evaluation is fast and reliable also near the logarithmic 
singularity at the threshold. 
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FIG. 2. Fit to the data from Ref. [12] with a parametrized 
Im/. In the inset the resulting Im/ is shown compared with 
ImFo averaged over 2.3 < Q < 3.2 and properly scaled. 

We take advantage of a symmetry in Eq. (Q) to define 
a dimensionless function f{uj) = XF(uj), constrained by 
the normalization condition dujlmf{uj) = ujn — 2A. 

We thus define gs = Vs/X^^^, and /3 = P/{\^^^a). In this 
way a, f3, gs, and the N — 1 parameters that define / 
are independent and can be fitted to the data. We set 
Im/(w) = qn for uj > lun and subtract out the infinite 
constant that would appear in Rei^(w). 

We thus fitted Eq. (convolved with the known in- 
strumental resolution) to experimental data of Ref. | |T^ ] 
(at 1.30 K) by minimizing numerically for tti^ param- 
eters: a, f3, 53, w(Qi), . . ., w(Q„), ai, . . ., un-i (N = 15 
for the fit presented). The minimization procedure has 
been performed with different standard routines and the 
result turned out to be independent of the choice of {tOi} 
and the initial value of a^. A typical starting point for 
Im/ is simply ai — 1 for all i. The resulting fit is shown 
in Fig. 1^. Good agreement between theory and experi- 
ment is obtained with a reduced (Xr) of nearly 4, thus 
indicating that even if we are leaving a large freedom in 
/ the agreement is significant. 

Our main result is summarized by the dispersion rela- 



tion for the undamped excitation shown in Fig. ^ where 
it is compared with the one reported in Ref. |1J] . We find 
that the model can quantitatively explain that the peak 
position of S{Q, uj) is slightly larger than 2A for Q > 2.6 
A~^. This originates from a mixing within the instru- 
mental resolution of the contribution of the sharp peak 
at energy luq slightly smaller than 2A with that of the 
continuum of two rotons excitations starting at 2 A. On 
general grounds the continuum should depend strongly 
on UJ near 2A, as for uj > 2A there are much more states 
available to decay into. For these reasons, the assump- 
tion of a background of constant slope, used to obtain the 
values of Ref. p^ , is not valid in this momentum region. 
Our procedure exploits instead the theoretical model for 
x{p) to extract the value of ujq. In this way we can find 
the final part of the dispersion relation with improved ac- 
curacy and it turns out that data agree with a dispersion 
relation for the excitations always below 2A. 
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FIG. 3. Dispersion relation found in the present work (solid 
line), and in Ref. [4,14] (dashed line). The inset shows the 
complete dispersion relation. 

The importance of the fitted parameter /? has been 
checked by studying the function x^(/5); where all the 
other parameters are properly modified to minimize 
for each value of /3. The confidence region for /3, i.e. 
the values of (3 such that X^{P)/Xmin < l-^: turns out 
to be -0.17 < f3 < -0.01 meV"^/^ with a best value 
(3 — —0.06 meV~^/^. This implies that the direct ex- 
citations of two quasiparticles by the neutron gives a 
small but sizable contribution to ^(Q,^). Concerning 
the other parameters we find = 1.4 and 173 ~ 0.8 
meVi/2. 

The shape of Im/ found by the fit has two main 
features: a clear peak at uj ~ 2 meV and a "quasi- 
divergent" behavior at the threshold. The peak is due 
to the maxon-roton van Hove singularity. This can be 
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verified by comparing the function lmFo{Q, ui), averaged 
over 2.3 < Q < 3.2 A~^and properly scaled, with the 
fitted Im/(w) (see the inset of Fig. ||). Hence it is clear 
that the peak corresponds in shape and position to the 
maxon-roton singularity. It is remarkable that no trace 
of the peak is apparent in any of the experimental plots. 
It is only by exploiting the correlation between plots with 
different Q that we have been able to extract this infor- 
mation. On this basis, using the fitted parameters it is 
also possible to predict that a peak and its shape should 
be observable with a resolution of 0.1 meV (to be com- 
pared with 0.5 meV of Fig. ||). 

The quasi-singularity at threshold can be understood 
as an interaction effect, namely a signature of the attrac- 
tive interaction between rotons. As a matter of fact, in 
the small region of energy near the threshold we can ap- 
ply Pitaevski-ZRS theory to evaluate F{uj) that 
reads: 

F(Q,c^)=4pln— ^ (^l-54ln— , (5) 

where p — — Imi^o(Q, ^ — 2A+)/27r is the threshold den- 
sity of states, 54 = 2V4P, and Z? is a cutoff that can 
be set to 1 meV as changes in D can be easily reab- 
sorbed in small variations of ujo, gs, and 54. One can 
thus readily verify that with a negative value of 94 Eq. (|^) 
gives a ImF that for lu 2A+ reproduces qualitatively 
Im/ of Fig. ^. To verify quantitatively this fact and 
to find an estimate of 54 we have repeated the data 
fit using Eq. (||) to parametrize /(w), (instead of the 
spline parametrization) and setting a cutoff in energy at 
2A-I-0.2 meV, in such a way to apply Eq. ^ only where 
it is supposed to hold. The resulting Xb, — 1-7 for the 
fit gives evidence that the theory works quantitatively 
in this region and we get for the interaction parameter 
= .g4/(2p) « -4.7 meV (54 = -0.15). The bound- 
state energy that we obtain from Eq. (||) is indeed very 
small Eb = D exp{l/g4} « 1.3fieV. This second fit 
gives also an additional estimate of the dispersion re- 
lation that agrees with the previous one and confirms 
that an undamped state is present up to Q — 3.6 A^^. 
The confidence region for 34 is —0.22 < 34 < —0.08, re- 
stricting P to its confidence region already found, which 
indicates that the interaction is definitely attractive. 

Some final comments are in order. First, while Fig. |^ 
shows the fit for 2.3 < Q < 3.2, the inclusion of the addi- 
tional set of data for Q — 3.6 increases slightly Xr (since 
the momentum range over which we are assuming the 
vertices to be momentum-independent may be too large) 
but it remains in any case a good fit to data. For this 
reason we report in Fig. |^ the value for cij(Q) obtained 
in this way. Second, we studied the cutoff dependence of 
x\ when /(w) is parametrized according to Eq. (||) and 
we found that it is very weak up to w = 2 meV where 
the effect of the maxon-roton peak becomes important. 
Thus use of the Pitaevskii-ZRS theory for F in our ex- 
pression (0) does not give a good description of the data 



if the cutoff in energy is removed. This indicates that the 
maxon-roton structure and, in general, the whole shape 
of Im/ plays a crucial role in determining ^(Q,^). 

In conclusion, we presented a theory for S'(Q,cj) that 
takes into account both one- and two-quasiparticle exci- 
tations by the neutron. The theory reproduces the exper- 
imental result over a large range of energy and momen- 
tum. We have thus been able to extract the final part of 
the spectrum dispersion relation in ^He. Our theory can 
be regarded as an extension of the Pitaevskii-ZRS theory 
taking into account the effect of two-particle excitations. 
Moreover the range of validity is enlarged as we make no 
hypothesis on the F function, but we extract it directly 
from data. In the region where Pitaevskii-ZRS theory 
holds we have used it to parametrize F in Eq. (Q) and we 
found that the interaction potential among rotons (V4) 
is attractive in this momentum region. 

I am indebted to P. Nozieres for suggesting this prob- 
lem to me and for many discussions. I gratefully acknowl- 
edge B. Fak for discussions and critical reading of the 
manuscript. I acknowledge B. Fak and J. Bossy for let- 
ting me use their data prior to publication. I also thank 
N. Cooper, N. Manini, A. Wiirger, P. Fieri, G.C. Strinati, 
and H.R. Clyde for useful discussions. 



* electronic address: pistoles@ill.fr. 

[1] A. GrifRn, Excitations in a Bose- Condensed Liquid 
(Cambridge University Press, Cambridge, UK, 1993). 

[2] H. R. Clyde, Excitations in Liquid and Solid Helium (Ox- 
ford University Press, Oxford, UK, 1994). 

[3] L. P. Pitaevskii, Sov. Phys.-JETP 36, 830 (1959). 

[4] R. A. Cowley and A. D. B. Woods, Can. J. Phys. 49, 177 
(1971). 

[5] J. Ruwalds and A. Zawadowski, Phys. Rev. Lett. 25, 333 

(1970). A. Zawadowski, J. Ruvalds, and J. Solana, Phys. 

Rev. B 5, 399 (1972). 
[6] L. P. Pitaevskii, JETP Lett. 12, 82 (1970). 
[7] K. Bedell, D. Pines, and A. Zawadowski, Phys. Rev. B 

29, 102 (1984). 
[8] K. J. Juge and A. GrifRn, J. Low. Temp. Phys. 97, 105 

(1994). 

[9] A. J. Smith et al, J. Phys. C 10, 543 (1977). 

[10] D. K. Lee, Phys. Rev. 162, 134 (1967). 

[11] B. Fak and K. H. Andersen, Phys. Lett. A 160, 469 
(1991). B. Fak, L. P. Regnault, and J. Bossy, J. Low 
Temp. Phys. 89, 345 (1992). 

[12] B. Fak and J. Bossy, J. Low Temp. Phys. in press (1998). 

[13] Within the standard notation for the Up and Hp func- 
tions we find in the Bogolubov approximation: a(p) = 
nV^iup - Vp), /3(p,q) = -«p+q«q, and 7(p,q) = 
Up+qUq -I- VdVp+d, with Uo the condensate density. 

[14] R. J. Donnelly, J. A. Donnelly, and R. N. Hills, J. Low 
Temp. Phys. 44, 471 (1981). 



4 



